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Do Thermoelectric Materials in Nanojunctions Display Material Property or Junction 

Property? 

Yu-Chang Chc^ and Yu-Shcn Liu 
Department of Electrophysics, National Chiao Tung University, 1001 Ta Hsueh Road, Hsinchu 30010, Taiwan 

The miniaturization of thermoelectric nanojunctions raises a fundamental question: do the ther- 
moelectric quantities of the bridging materials in nanojunctions remain to display material properties 
or show junction properties? In order to answer this question, we investigate the Seebeck coefflcient 
S and the thermoelectric figure of merit ZT especially in relation to the length characteristics of 
the junctions from the first-principles approaches. For S, the metallic atomic chains reveal strong 
length characteristics related to strong hybridization in the electronic structures between the atoms 
and electrodes, while the insulating molecular wires display strong material properties due to the 
cancelation of exponential scalings in the DOSs. For ZT, the atomic wires remain to show strong 
junction properties. However, the le ngth ch rematistics of the insulation molecular wires depend on 
a characteristic temperature To = \/W/j(J) around 10 K. When T <C To, where the electron trans- 
port dominates the thermal current, the molecular junctions remain to show material properties. 
When T 3> To, where the phonon transport dominates the thermal current, the molecular junctions 
display junction properties. 
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Nanoscale thermoelectric devices can be considered as 
the new types of devices which can be embedded into 
integrated chip sets to assist the stability of devices 
by converting the accumulated waste heat into useful 
electric energy. There has been an ever increasing in- 
terest in the thermoelectric properties of nanojunctions 
[II, i, S, i, H, S 0, i, i, [3, partially motivated by 
the recent experiments demonstrating the capability of 
measuring the Seebeck coefficients in the molecular junc- 
tions [111, As the Seebeck coefficients are relevant 
not only to the magnitude but also to the slope of den- 
sity of states (DOSs), they can reveal more detailed in- 
formation about the electronic structures of the mate- 
rials sandwiched in the nanojunctions beyond what the 
conductance measurements can provide. The Seebeck 
coefficients have been applied to explore the electronic 
structures of molecular junctions using functional sub- 
stitutions for the bridging molecules [l( 
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have proposed using gate fields and external biases as 
means to modulate the Seebeck coefficients in nanojunc- 
tions 0, [3, [13] ■ Much research has been devoted to the 
study of Seebeck coefficients fl, 2, 3, 4, 1, [3, [3, S, [3, [3 , 
but little is known about the fundamental thermoelectric 
properties in nanojunctions. 

Indeed, the miniaturization of thermoelectric nano- 
junctions raises a fundamental question: does bridging 
thermoelectric material in nanojunctions show material 
properties or junction properties? Thermoelectric bulk 
crystals usually show material properties, where the ther- 
moelectric physical quantities arc irrelevant to the sizes 
and shapes of materials. In addition, recent experiments 
on Seebeck coefficients in molecular junctions also re- 
veal strong signals of material properties. These exper- 
iments have observed that Seebeck coefficients are in- 
sensitive to the number of molecules in junctions and 
show rather weak dependence on the lengths of the bridg- 
ing molecules, which is in sharp contrast to the conduc- 
tance which shows str ong exponential dependence on the 
lengths of molecules [ll[, [13, [13 ■ However, the bridg- 



ing materials in nanojunctions may have strong inter- 
actions with the contacts. From this point of view one 
can say thermoelectric quantities can display the junc- 
tion characteristics. Considering the examples and the 
reason quoted above, it is therefore not obvious whether 
the thermoelectric quantities of the bridging materials 
in nanojunctions display junction properties or material 
properties. 

In this letter, we will show that the thermoelectric 
quantities in nanojunctions unnecessarily display entire 
material properties or junction properties. To demon- 
strate this point, this study investigates two important 
thermoelectric quantities, the Seebeck coefficient (5) and 
the thermoelectric figure of merit (ZT), in (metallic) alu- 
minum atomic junctions and (insulating) molecular junc- 
tions. It shows that metallic atomic chains reveal strong 
junction properties while the insulating molecular wires 
partially possess the material properties, where S reveals 
the material property and ZT displays the junction prop- 
erty at temperatures larger than the characteristic tem- 
perature To- 

To answer the question, we have developed a theory 
with analytical expressions for S and ZT allied to a fully 
self-consistent first-principles calculation in the frame- 
work of the density functional theory (DFT). It allows 
us to numerically calculate S and ZT and subsequently 
investigate these quantities analytically. We focus on the 
subject on whether S and ZT depend on the character- 
istics of the junctions, especially on the dependence on 
length-characteristic of junctions. Before turning to the 
detailed discussion, let us begin with a brief introduction 
on how to calculate S and ZT. First, we consider that 
the junction consists of source-drain electrodes, with dis- 
tinct chemical potentials Hl{r) temperatures Ti^[r), 
as independent electron and phonon reservoirs. When an 
additional infinitesimal temperature AT is applied across 
the junction, an extra voltage AV is induced to com- 
pensate the electric current induced by the temperature 
gradient AT across the junction. We then derive the 



2 



expressions for S (defined as S ~ AV/AT) and the elec- 
tron thermal conduetance (defined as Kei = AJq / AT), 
where Jg is the thermal current conveyed by the elec- 
trons which also carry the electric current). 
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the transmission function t{E) ~ t^{E) = t^{E), 
which is a direct consequence of the time-reversal 
symmetry. It has been assumed that the left and 
right electrodes serve as independent electron and 
phonon reservoirs where the electron population is 
described by the Fermi-Dirac distribution function, 

/^(^^ = 1/ [cxp {{E - fiLiR)) IksTLiR)) + 1] , and ke 
is the Boltzmann constant. The transmission functions 
are computed using the wave functions_obtained self- 
consistently in the DFT framework IJ, 
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notice that the above equations are suitable for describing 
S and in nanojunctions operated under finite external 
biases, where two electrodes can have different tempera- 
tures. 

In addition, the differential conductivity, typically in- 
sensitive to temperature in cases where direct tunneling 
is the major transport mechanism, can be expressed as 



i=L,R 



- fE)riE)/kBT,. (3) 



So far, the physical quantities {S, a, and k^i) which 
have been discussed are related to the electron transport 
. It must be noted that the heat current is conveyed by 
the electrons and phonons simultaneously. The phonon 
thermal conductance (Kph) usually dominates the com- 
bined thermal conductance at large tem- 
peratures. The complete discussion on ZT shall include 
the essential ingredient Kph] thus ZT can be expressed 
as follows: 



ZT = 



Kel + Uph 



(4) 



where T = {Tl + Tb)/2 is the average temperature of 
the source-drain electrodes. We estimate the phonon 
thermal conductance following the approaches of Patthon 
and Geller [l^. It is assumed that the nanojunction is a 
weak clastic link, with a given stiffness which can be eval- 
uated from the total energy calculations, attached to the 
electrodes modeled as phonon reservoirs. The phonon 
thermal conductance (defined by kpt — AJq'/AT) is 
given by: 
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where n^ji^ = l/[cxp{E/kBTL,^ji^) - 1] and NL(^ii){E) « 
C • i? is the Bose-Einstein distribution function and the 
spectral density of phonon states in the left (right) elec- 
trode, respectively. The stiffness of the bridging nanos- 
tructure is = YA/l, where Y is the Young's modulus 
and A {I) is its cross-section (length). 

We have numerically computed S and ZT using 
Eqs. H]) to ([5]) allied to the transmission functions ob- 
tained self-consistently in the DFT framework, as shown 
in Figs. [T] and [51 To elaborate the properties of the See- 
beck coefficient and the the thermoelectric figure of merit, 
we will limit our discussion to the linear response regime 
(i.e., fiL ~ fJ-R = fJ-) and T^ = Tr = T. After expanding 
5, Keh E^nd Kph in terms of the temperature T, we obtain 
the analytical expressions for the S and ZT. 



S ~ aT, 



(6) 



where a = — Tr^fc^ ^g^'' / (3er(/i)). We have noted that 
the Seebeck coefficient depends on the magnitude and 
the slope of the transmission function, and is linearly 
proportional to T at low temperatures. 
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where we have expanded k^i and Kph up to the lowest or- 
der in temperatures as k^i ~ f3T and Kph ~ ^{l)T^ . The 
prefactor j3 and ^{l) are /3 = 2n^kgT{fi)/ {3h) and 7(/) = 
S-K^ k%C^ A'^Y^ / {15hP) , respectively. One may notice 
that there is a characteristic temperature To = 13/^(1), 
which is around 10 K for the alkanethiol molecular junc- 
tions and is negligibly small for aluminum atomic junc- 
tions. When T ^ To, the electron thermal conduc- 
tance dominates and ZT « aS'^T/kei « [a^cr//?] T^, 
which is irrelevant to the length-characteristic of the 
junction and is proportional to T^ as temperatures in- 
crease. When T ^ To, the phonon thermal conductance 
dominates and ZT tends to have a saturation value of 
ZT w aS^T/Kph « c?al^{J), which is related to the 
length of the junction. The above analytic expressions 
provide a convenient means for analyzing the length char- 
acteristic of S and ZT in nanojunctions. 

Before turning to the detailed discussion on the ther- 
moelectric properties in nanojunctions, we must draw at- 
tention to the reason why nanojunctions display material 
properties for S and ZT . The Seebeck coefficients are di- 
rectly related to the transmission functions, which are de- 
termined by the electronic structures of the crystal mate- 
rials irrelevant to the size and shape of the material, lead- 
ing to the material properties for S. Furthermore, the 
conductance u and the combined thermal conductance 
K = Kei + ^ph are proportional to the contact surface 
and inversely proportional to the length scale in the bulk 
crystal materials, which leads to the material properties 
for ZT due to the cancelation of the geometric factors 
in the conductance a and the combined thermal conduc- 
tance K. However, whether the thermoelectric materials 
in nanojunctions remain to display material properties or 
show junction properties have not yet been fully realized. 



3 




(c) 


^° T,K, 8° 


1.2x10^ 










c, 


w 8.0x10' 













& 4.0x10 ' 












0.0 







120 



10' 







C, 












c_ 




10 100 " 




T(K) 



















120 



10 

T(K) 



100 



FIG. 1: (color online) Alkanethiol junctions at Vs = 0.01 V 
for C4 [solid (black) lines] , Ce [dotted-dash (red) lines] , and Cs 
[dash (blue) lines]: (a) the Seebeck coefficient S vs. T; (b) the 
electron thermal conductivity n^i and the electric conductance 
(T (inset) vs. T; (c) the phonon thermal conductance Kp/i vs. 
T; and (d) the log-log plot of ZT vs. T (inset shows the case 
of neglecting the thermal conductance = 0) . 
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FIG. 2: (color online) Aluminum atomic junctions at Vb ~ 
0.01 V for 1 Al [solid (black) lines], 2 Al [dotted-dash (red) 
lines], 3 Al [dash (blue) lines], and 4 Al [dot-dot-dash (green) 
lines]: (a) the Seebeck coefficient S vs. T and the schematic 
of 3-Al atomic chain and its energy diagram (inset) where the 
Al — Al bond distance is about 6.3 a.u.; (b) the electron ther- 
mal conductance k^i and the electrical conductivity a (inset) 
vs. T; (c) the phonon thermal conductivity Uph vs. T; (d) the 
log-log plot of ZT vs. T (inset shows the case of neglecting 
the thermal conductance ~ 0). 



Now, let us look closely into the properties of S and ZT 
in the Unear response regime at low temperatures with 
T = Tl = Tr for the two catalogs of nano junctions: 
the (insulating) alkanethiol junctions and the (metallic) 
aluminum atomic junctions. 

First; we examine the dependence of S and ZT on the 
lengths of the alkanethiol junctions Alkanethiols 
[CH3(CH2)„-iSH, denoted as C„] are a good example of 
reproducible junctions which can be fabricated [la;[l§]. 
It has been well established that the non-resonant tun- 
neling is the main conduction mechanism in the C^i junc- 



tions. As shown in the inset of Fig. [H^b), the con- 
ductance is small, insensitive to temperatures, and de- 
creases exponentially with the length of the molecules, as 
a = Co exp (—££), where I is the length of C„ molecule and 
^ « 0.78 A-i [2fl,[2l|,[23,[2i,[2i]. By exploiting the peri- 
odicity in the {CH2)2 group of the C„ chains, the wave 
functions of the C„ junctions are calculated by a simple 
scaling argument, leading to exponential scaling in the 
transmission functions t{E). It must be noted that the 
Seebeck coefficients are proportional to Con- 
sequently, the Seebeck coefficient does not show depen- 
dence on the length [see Fig. [Ija)] of the C„ junctions 
due to the cancelation of the exponential scaling behav- 
ior in r(/i) for insulating molecules. We have noted that 
the cancelation may not be complete in the real experi- 
ments; thus, the Seebeck coefficients could possibly show 
weak length dependence due to other effects [ll|, ■ In 
such case, the Seebeck coefficients mostly display mate- 
rial properties. 

Let us now investigate the length-dependence of ZT 
for the C„ junctions. We predict that the C„ junctions 
display the material properties when T <^ Tq and display 
the junction properties when T ^ Tq. The reason for this 
is due to different scaling behaviors on the length char- 
acteristics of the junctions for the electron and phonon 
thermal conductance and the competition between them. 
When T <C To, the electron transport dominates the 
thermal current, and thus ZT « aS^T/ Kei- As shown 
in the inset of Fig. [Hd), ZT is independent from the 
lengths of C„ molecules and displays material proper- 
ties. This result can be explained quite naturally by the 
cancelation of the exponential scaling in a and k^i be- 
cause both the electric current (with conductance a) and 
the electron thermal current (with thermal conductance 
Ke/) are conveyed by electron transport. When T ^ Tq, 
the phonon transport dominates the thermal current and 
thus ZT « aS'^T/ Kpii- As shown in the main body of 
Fig- Hd), ZT oc Z^exp(— ^Z) because of cr oc exp(— ^Z) 
[see the inset of Fig.[TJb)], and Kph oc [see Fig.[T][c)]. 
In this case, ZT displays junction properties. 

Finally, we will examine the dependence of S and ZT 
on the lengths of the metallic aluminum (Al) atomic 
jimctions [l7|. An Al atomic chain is an ideal testbed 
for studying the charge transport at the atom-scale level 
[see the inset of Fig.J^a) for a schematic of the aluminum 
atomic junction] This study has observed 

that S and ZT depend on the geometric characteristic of 
the junctions and show the junction properties in any 
case, which is in sharp contrast to the C„ junctions. At 
a fixed temperature, it has been observed that ZT and 
the magnitude of S increase as the number of Al atoms 
increases, as shown in Figs. [2ja) and (d). The increase 
of the Seebeck coefficients is due to the increase of the 
slopes in the transmission functions at the Fermi levels. 
The reason for this may be due to the strong hybridiza- 
tion between the electronic structures of atoms and the 
electrodes. The negative sign of the Seebeck coefficients 
indicates that the metallic Al atomic junctions arc n- 
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type [the Fermi energy is closer to the lowest unoccupied 
molecular orbital (LUMO)]. 

In conclusion, this study has raised and answered an 
important fundamental question in thermoelectric nano- 
junctions: do the thermoelectric materials in nanojunc- 
tions remain to display material properties or show junc- 
tion properties? To answer this question, we have devel- 
oped a theory with analytical expressions for S and ZT 
allied to a fully self-consistent first-principles calculation 
in the DFT framework. Using the insulating alkanethiol 
molecular junctions and metallic aluminum junction as 
examples, this study concludes that the thermoelectric 
quantities in nanojunctions do not necessarily display 
material properties or junction properties. The metallic 
atomic chains reveal strong length characteristics related 
to the strong hybridization in the electronic structures 
between the atoms and electrodes, while the insulating 
molecular wires display strong material properties due to 
the cancelation of exponential scalings in the DOSs. For 
S, the metallic atomic chains reveal strong length char- 
acteristics related to the strong interactions between the 
electronic structure of atoms and the electrodes, while 
the insulating molecular wires show independence from 
the lengths of molecules, thus displaying strong material 
properties due to the cancelation of exponential scalings 
in the density of states. It may be worth pointing out 
the cancelation may not be complete in real experiments 
due to other effects; the important point is, the strong 
exponential scaling behavior with the lengths of the junc- 
tions should be canceled. For ZT, the atomic wires re- 
main to show strong junction properties. However, the 
length characteristics of the insulation molecular wires 
depends on a characteristic temperature Tq — \/ (3/^{l) 
around 10 K. When T Tq, where the electron trans- 
port dominates the thermal current, the molecular junc- 
tions remain to show material properties. When T ^ Tq, 
where the phonon transport dominates the thermal cur- 
rent, the molecular junctions display junction properties. 
The length characteristic of the C„ molecules is according 
to ZT (X P exp (— CO- The different length-scaling behav- 
iors between the thermal current conveyed by the elec- 
trons and the thermal current conveyed by the phonons 
offer the key to the understanding of characteristic tem- 
perature and the behavior of ZT. We believe that this 
study is a substantial step towards the understanding of 
the thermoelectric properties in nanojunctions, and we 
hope that this study will generate more experimental and 
theoretical explorations in the properties of thermoelec- 
tric nanojunctions. 
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